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We analytically investigate transverse magnetic (TM) spatial bright solitons, as exact solutions of 
Maxwell's equations, propagating through nonlinear metamaterials whose linear dielectric permit- 
tivity is very close to zero and whose effective nonlinear Kerr parameters can be tailored to achieve 
values not available in standard materials. Exploiting the fact that, in the considered medium, linear 
and nonlinear polarization can be comparable at feasible and realistic optical intensities, we identify 
two novel self-trapping mechanisms able to support two-peaked and flat-top solitons, respectively. 
Specifically, these two novel mechanisms are based on the occurrence of critical points at which the 
effective nonlinear permittivity vanishes, the two mechanisms differing in the way the compensa- 
tion between linear and nonlinear polarization is achieved through the non-standard values of the 
nonlinear parameters. 

PACS numbers: 78.67.Pt, 42.65.Tg 



I. INTRODUCTION 

Metamaterials have recently attracted a large atten- 
tion since they are characterized by unconventional elec- 
tromagnetic properties and have recently suggested novel 
and fascinating applications [IHEl • I n the context of non- 
linear optics, metamaterials have suggested two different 
ways to overcome its fundamental limit consisting in the 
fact that, due to the weakness of the nonlinear response, 
nonlinear behavior is observed at large optical intensi- 
ties. In the first scheme the material is tailored in such 
a way that the nonlinear medium is placed at the region 
where a large field enhancement occurs @, 0], whereas 
in the second scheme the linear dielectric permittivity is 
chosen to be so small to allow the nonlinear polariza- 
tion to fully rule the electromagnetic behavior {extreme 
nonlinear regime) @. Such an extreme nonlinear regime 
has been recently allowed to predict genuinely novel and 
intriguing phenomena as, for example, beam transverse 
power flow reversing Q where Poynting vector flips its 
sign along the transverse profile of an optical beam since 
optical radiation nonlinearly produces a region (around 
the beam propagation axis) where the dielectric behaves 
as a metal. A nonlinear slab characterized by a very small 
optical intensity has been shown to support transmis- 
sivity directional hysteresis [loj and the corresponding 
underlying electromagnetic multiplicity has been shown 
to be produced by the nonlinear field matching at the 
slab interfaces in the extreme nonlinear regime. In the 
context of nondiffracting and self-trapped waves, nonlin- 
ear epsilon-near-zero metamaterials have been recently 
shown to be able to support dark solitons characterized 
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by a central and extended region where Poynting vec- 
tor exactly vanishes together with the transverse electric 
field component [ill ]. 

In this paper, we analytically investigate trans- 
verse magnetic (TM) spatial bright solitons propagat- 
ing through epsilon-near-zero nonlinear Kerr metama- 
terials (material proposed in Ref.Q) and we predict 
the existence of two novel perfect soliton [l2[ families 
whose Poynting vector profiles are two-peaked and flat- 
top [13, [lj] , respectively. Such metamaterials host the 
extreme nonlinear regime where, due to the small value 
of the linear dielectric permittivity, linear and nonlin- 
ear polarizations can be comparable at feasible optical 
intensities. As a consequence, critical electromagnetic 
points are allowed where the overall effective nonlinear 
permittivity vanishes, and in this paper we prove that 
these points can be profitably exploited to tailor the soli- 
ton shape, allowing us to obtain both two-peaked and 
flat-top solitons. The occurrence of such points literally 
provides novel physical self-trapping mechanisms (much 
more involved than standard compensation between self- 
focusing and diffraction) where the longitudinal electric 
field component (i.e. the component parallel to the soli- 
ton propagation direction) is the basic mean to produce 
the above mentioned balance between linear and nonlin- 
ear polarizations. The full exploitation of such novel self- 
trapping mechanisms is possible since the nonlinear pa- 
rameters of the metamaterial Kerr response (as discussed 
in Ref. [8] ) can be tailored to assume values not available 
in standard materials. In order to prove the feasibility 
of the two novel soliton families, we propose a realistic 
composite medium characterized by a sufficient number 
of parameters which can be independently tuned to al- 
low the overall homogenized medium response to meet 
the conditions required for observing both two-peaked 
and flat-top solitons. In addition, as a consequence of 
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the small value of the linear permittivity, two-peaked and 
flat-top solitons turn out to be observable at feasible op- 
tical intensities. 

The paper is organized as follow. In Section II we 
discuss epsilon-near-zero nonlinear metamaterials, we in- 
troduce the basic theory for investigating bright solitons 
and we derive their existence conditions even for values 
of the nonlinear parameters not available in nature but 
available through the proposed nonlinear metamaterial. 
In Section III and Section IV we consider two-peaked 
and flat-top solitons, respectively, and we focus on their 
own underlying self-trapping supporting mechanisms by 
discussing their mutual differences and the role played 
by the unconventional values of the nonlinear parame- 
ters. In Section V we propose a feasible and realistic 
composite medium (nonlinear metamaterial) suitable for 
observing two-peaked and flat-top solitons. Finally, in 
Sec. VI, we conclude with a discussion and summary of 
our results. 



II. BRIGHT SOLITONS IN NONLINEAR 
EPSILON-NEAR-ZERO METAMATERIALS 

Consider monochromatic transverse magnetic (TM) 
fields E = [Eexp (-iuit)], B = Re [B exp (—iwt)] (of fre- 
quency uj), where 



E = E x (x, z)e x + E z {x, z)e z 
H = H y (x,z)e y , 



(1) 



propagating through a nonlinear metamaterial whose ef- 
fective medium response is described by the constitutive 
relations 

D = e eE + e oX [|E| 2 E + 7 (E • E)E*] , 

B = a^omH, (2) 

where e and fj, are the linear dielectric permittivity and 
magnetic permeability, respectively, whereas \ an d 7 are 
parameters characterizing the overall effective nonlinear 
Kerr response. As opposed to conventional media where 
e > 1, /i ~ 1 and 7 assumes three values depending on 
the physical mechanism supporting the Kerr nonlinear 
response [l5[ (i.e., 7 = for electrostriction, 7 = 0.5 for 
nonresonant electronic response and 7 = 3 for molecular 
orientation), we hereafter focus on suitable metamate- 
rials designed to exhibit a very small dielectric permit- 
tivity |e| <C 1 [16], an arbitrary \i and with 7 attain- 
ing, in principle, any real value. The possibility of tai- 
loring both linear and nonlinear medium responses has 
been investigated in the pioneering paper by J. E. Sipe 
et al. of Ref. [17( where the authors consider a composite 
structure comprised of spherical inclusion particles em- 
bedded in a host material. Remarkably, if the inclusions 
respond linearly and the host material responds nonlin- 
early, the composite is characterized by a nonlinear ef- 
fective Kerr response whose parameter 7 can span the 
range < 7 < 3 [13 ■ Recently, A. Ciattoni et al. [|| 



have improved these results by showing that TM optical 
waves travelling through a nonlinear layered medium, in 
the long wavelength limit, experience a nonlinear opti- 
cal response described by Eqs.© whose electromagnetic 
parameters can be independently tailored. In particular, 
using positive and negative (linear and nonlinear) dielec- 
tric layers, they have shown that a suitable composite 
tailoring allows 7 to achieve any prescribed value (even 
encompassing negative values) and e, at the same time, 
to assume a value very close to zero. The evident con- 
sequence of operating with a very small dielectric per- 
mittivity is that the extreme nonlinear regime can be 
observed where the nonlinear polarization can be com- 
parable with the linear part Q, i.e. the condition 



x[|E| 2 E + 7 (E-E)E*] «eE 



(3) 



can be satisfied at small and feasible optical intensities. 
Note that in standard media where e w 1, condition of 
Eq.Q is strictly meaningless since the Kerr nonlinear- 
ity is always a perturbation to the linear polarization 
and, for very high intensities, it always displays a form 
of saturation thus departing from the cubic behavior and 
never reaching values of the order of eE. Therefore, 
even though Eqs.([2]) formally coincide with the standard 
nonlinear Kerr response, the considered media host a 
marked and hitherto unexplored nonlinear behavior (ex- 
treme nonlinear regime) accessible at feasible optical in- 
tensities. It is worth stressing that, since the very small 
value of the permittivity is here achieved by averaging 
negative and positive dielectric permittivities, it is evi- 
dent that suitable gain media have to be inserted among 
the underlying medium constituents in order to compen- 
sate losses due to the presence of negative dielectrics (see 
Sec. V). 

We focus on non-diffracting solitary waves propagating 
along the z-axis of the kind 



(4) 



where £ = ^\[ie\k x, £ — y/\fie\koZ (dimensionless spa- 
tial coordinates), fco = oj/c and we require both j3 (di- 
mensionless propagation constant) and u x and u z (di- 
mensionless field amplitudes) to be real quantities. Note 
that, using dimensionless field amplitudes introduced 
through Eqs.Q, the extreme nonlinear regime condition 
of Eq. (J3j) can be stated as 
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1. 



(5) 



Maxwell equations V x E = io;B, V x H = — iusT) for the 
fields of Eqs.([T|) and Eqs.Q and with the constitutive 
relations of Eqs.([2]) yield, after some algebra, the system 
of first order differential equations 



0- 



, du z 



/3 2 



(NL) 



du x 



(Q\ TABLE I: Bright soliton existence ranges of u^. depending 



where a e — ex/l e xl; a n — A'X/ImxI an d we have defined 

ef^ = a e + (l+ 7 )^ + (l-7)^, 

4^) = a £ + (l- 7 )^ + (i + 7 ) M 2 j (7) 



as effective normalized nonlinear dielectric permittivi- 
ties. We note that the system of Eqs.© is equivalent 
to Maxwell's equations if the condition 



+ u x f 0, 



(8) 



holds along the whole field profile. It is worth stressing 
that the system of Eqs.© is integrable since it admits 
the first integral 

F(u Xl U z ) = [(ft 2 - CTeV^ul - CTeO-^ul} 



-^{(/3 2 - a £ a„) - <r M [(l + j)u 2 x + (1 - 7 )«»]} 



2 2 



(9) 



i.e. the relation dF/d£ — is satisfied along the solu- 
tions u x (^), u z (£) of Eqs.© (see Ref.Q for details). We 
here investigate bright solitons whose normalized electric 
field components u x and u z asymptotically vanish at in- 
finity and are spatially even (u x (£) = u x (— £)) and odd 
(u z (0 = ~ w z(~0): so that we consider the boundary 
conditions u x (0) = u x0 , u z (0) = 0, u x (+oo) — and 
u z (+oo) = 0. The vanishing of the field at infinity im- 
plies that the first integral F has to vanish along the 
whole soliton profile, i.e. F(u x (£), u z (0) = 0, a rela- 
tion which, evaluated at £ = 0, allows us to analytically 
evaluate the propagation constant 
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2(7, 



[a e + (1 + tKq 
2a e + 3(l+j)ul 



(10) 



Remarkably, bright solitons are represented in the phase 
space (u x ,u z ) by the special level curves F{u x ,u z ) = 
which are orbits of the dynamical system of Eqs.© 
(known as homoclinic orbits) joining the equilibrium 
point (u x ,u z ) — (0,0) to itself. The main consequence 
of this identification is that bright solitons asymptoti- 
cally coincide for £ — > ±oo with the linear inhomoge- 
neous transverse plane waves e^V '^ 2 - CT e°>£ ^ respectively, 
so that the relation 
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a £ a^ > 



(11) 



is a necessary condition for soliton existence. The exis- 
tence of at least one turning point (i.e. where du z /d£, = 



on 7, a E , cr M . The parameter a is defined in Eg. (fill) 



occurs) along the orbit corresponding to a bright soli- 
ton is an additional condition to be met and it is ful- 
filled by requiring the curves F(u x ,u z ) — and [f3 2 — 
&ti£x NL,> }u x = (see the first of Eqs.©) to have an inter- 
section. These necessary conditions for the existence of 
bright solitons, together with Eq.©, allow us to obtain 
the results of Table I where we report, for each possible 
combinations of a e and , the range of u 2 (depending 
on 7) where bright solitons can be excited. 



(a) 
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FIG. 1: Profiles of u x (£) (panel (a)), u z (£) (panel (b)) and 
normalized Poynting vector S z (£)/So (panel (c)) of various 
bright solitons with cr e = 1, <r M = 1, sign(x) = 1, = 2 and 
7 spanning the range —0.8 < 7 < 10. 



III. TWO-PEAKED SOLITONS 

A first group of bright solitons (see Table I) exists for 
any values u x0 in the range 7 > — 1 and for er c = 1 , = 1 
(i.e. for right-handed metamaterial with focusing non- 
linearity, % > 0, or for left-handed metamaterials with 
defocusing nonlinearity, \ < 0). Note that this group 
of bright soliton encompasses the perfect optical soli- 
tons described in Ref.[l2| where authors solely consid- 
ered the standard Kerr nonlinearity arising from the elec- 
tronic nonresonant response (7 = 0.5). In panels (a) 
and (b) of Fig.l we report the bright soliton profiles of 
u x and u z , respectively, for u x0 = 2, —0.8 < 7 < 10 
and cr e = l,o> = 1. For completeness, in Fig. 1(c) 
we report the profiles of the normalized Poynting vec- 
tor S z ( £)/S = sign(n)[ (f3u x - du z /d^)u x ]e z (where 
So = \J (e |e| 3 )/(4^olMlx 2 )) corresponding to the solitons 
reported in Fig. 1(a) and Fig. 1(b). From Fig.l, we note 
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FIG. 2: (Color online) Plots of F(u x , u z ) evaluated for <r e = 1, 
o> = 1, u x0 = 2, 7 = 0.5 (panel (a)), 7 = 3 (panel (b)). Black 
lines correspond to the homoclinic orbits whereas the stars 
and the circles are located at the saddle points and maxima 
of F(u x ,u z ), respectively. 



that two-peaked bright solitons can be excited, i.e. soli- 
tons whose electric field x-components and Poynting vec- 
tors show two pronounced peaks along their transverse 
profiles, the higher the value of 7 the more pronounced 
the two peaks. In order to grasp the underlying mech- 
anism producing such an atypical (for the Kerr nonlin- 
earity) two-peaked soliton profile, we consider the equi- 
librium points of the system of Eqs.((6]) (i.e. du x /d£ = 
and du z /dt; = and hence necessarily coinciding with 
the critical points of the first integral F(u x ,u z )) which, 
in addition to the origin (u x ,u z ) — (0,0), are given by 
the relations 



0. 



(12) 



Note that the number of the critical points of F depends, 
in addition to u x0 , on the value of the nonlinear param- 
eter 7. In fact, for I7I < 1 (see Fig. 2(a)) F has only a 
maximum at the point (u x ,u z ) — (y/ \(3 2 — 1)/(1 + 7), 0) 
(belonging to the axis u z = 0) and a saddle point at the 
origin (u Xl u z ) — (0,0); on the other hand, in the case 
7 > 1 and u x0 > S (see Fig. 2(b)) where 



£(7) 



(7 + 2) + ^5 7 2 + 2 7 
2( 7 2 - 1) 



(13) 



F has two saddle points (located at (u x ,u z ) = 
(V(/3 2 - !)/(! + 7), 0) and the origin (u x ,u z ) = (0,0)) 
and two maxima not belonging to the axis u z = (and 

hence characterized by the condition ei NL ^ = 0, see 
Eqs. (fT2"|) ). From Fig. 2 it is evident that, in the consid- 
ered situations, the homoclinic (solitonic) orbits make a 
turn around the maxima of the function F(u x , u z ) so that 
the shape of the orbits (and consequently the soliton pro- 
files) depends on the number of maxima enclosed by the 



homoclinic loop. Therefore, since there are one and two 
maxima for I7I < 1 and 7 > 1 (for u x0 > S), respectively, 
we conclude that homoclinic orbits undergo a qualitative 
structural change at 7 = 1. Correspondingly, for I7I < 1 
homoclinic orbits are simple loops (see Fig. 2(a)) and soli- 
tons are bell-shaped whereas for 7 > 1 (for u x0 > S) 
homoclinic orbits are structured loops (see Fig. 2(b)) and 
two-peaked solitons occur. In other words (in a language 
more close to nonlinear optics) the discussed two kinds 
of bright solitons are literally supported by two different 
self-trapping mechanisms, and the nonlinear metamateri- 
als employed in this paper provide the additional param- 
eter 7 to switch from one mechanism to the other. It is 
worth noting that the first of Eqs. (fl~2"|) globally holds for 
all the considered maxima whereas the second of Eqs.([T2|) 
is differently fulfilled since, for I7I < 1 and 7 > 1 (for 

u xQ > <5), maxima occur for u z = and ei NL ^ = 0, 
respectively. Therefore, bell-shaped soliton are not con- 
strained by additional specific requirements and they ac- 
tually coincide with perfect solitons of Ref. [l2[ . On the 
other hand, two-peaked solitons supported by the novel 
self-trapping mechanism are due (for u x0 > 5) to the 
existence of off-axis maxima of F corresponding to the 
condition e{ NL ^ = 0. However, the normalized effective 
nonlinear permittivity €z can vanish only if the non- 
linear contribution to the polarization (1— 7)it 2 .+(l+7)u 2 
exactly balances the linear part a e — 1 (see the second 
of Eqs.((7j)) and it is worth repeating that this condition 
(equivalent to Eqs.([3])) is realistically achievable solely 
by means of the nonlinear metamaterials with dielectric 
permittivity very close to zero we are considering in this 
paper. 



IV. FLAT-TOP SOLITONS 



A second group of bright solitons (see Table I) exists 
for u x0 < a where 



01(7) 



2(|7|-l)( V / l7l + l) 



(14) 



in the range 7 < — 1 and for cr e = — 1 and <r M = — 1 
(i.e for right-handed metamaterials with defocusing non- 
linearity, x < 0, or for left-handed metamaterial with 
focusing nonlinearity, \ > 0). In Fig. 3(a) and Fig. 3(b) 
we report the soliton profiles (for various u x0 ) of u x and 
u z , respectively and in Fig. 3(c) we plot the correspond- 
ing normalized Poynting vector S z /So, for sign{e) = 1, 
sign(fi) — 1, sign(x) = —1, and 7 = —1.5. Note that, 
for each 7 < —1, if u x0 is very close to a (i.e. for 
u xQ > 0.99a), flat-top bright solitons occur since their 
electric field x-component and Poynting vector profiles 
are characterized by a flat and sharp core region (where 
the intensity is almost constant) which abruptly stops 
being surrounded by lateral regions where the field ap- 
proximately vanishes. On the other hand, for smaller in- 
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FIG. 3: Profiles of u x (£,) (panel (a)), u z (£) (panel (b)) and 
normalized Poynting vector S z (^)/So (panel (c)) of various 



(with different 



bright solitons for <r e 



-1. cr u 



sign(x) 



-1 and 7 = —1.5. 



tensities (i.e. for < 0.99a) solitons are characterized 
by a bell-shaped profile. 




FIG. 4: (Color online) Plots of F(u x ,u z ) evaluated for a t — 
— 1, <t m = — 1, 7 = —7, u x0 = 0.8a (panel (a)), u x0 = 0.99a 
(panel (b)). Black lines correspond to the homoclinic orbits 
whereas the stars and the circles are located at the saddle 
points and maxima of F(u x ,u z ), respectively. 



Even the mechanism leading to the formation of flat- 
top solitons can be understood by considering the critical 
point of the first integral F. For 7 < — 1, it is simple to 
prove that the two off-axis critical points (i.e. with u z 7^ 

and satisfying t z NL ^ — 0) obtained by Eqs.(fl2]l are saddle 
points of F and this is particularly evident from Fig. 4(a) 
and Fig. 4(b) where we plot such saddle points (marked 
with the stars) together with the first integral F on the 
phase space (u x ,u z ) for u 2 = 0.8a and u 2 — 0.99a, 
respectively (for the parameters a e = —1, 07, = —1, 
7 = —7). Analogously to the case of the above discussed 
two-peaked solitons, the off-axis critical points of F have 
an incisive impact on the structure of the homoclinic or- 
bits and, in the case 7 < — 1 we are considering, the fact 



that they are saddle points implies that homoclinic orbit 
can not make a turn around them. Therefore, by varying 
u^o (for each 7 < — 1) the mutual position (and distance) 
between the homoclinic orbit and the saddle points can 
be changed (see Fig. 4) to the point that, for u 2 very close 
to a, the homoclinic orbit almost collides with the saddle 
points. When this happens, at the points (u x , u z ) of the 
homoclinic orbit which are very close to the saddle points, 
the RHS of Eqs.© are very small (since at the saddle 
they exactly vanish) and therefore du x /d^ and du z /d£, 
are very small as well, thus explaining the obtained flat- 
top soliton behavior. In the language of nonlinear optics, 
we can state that, for 7 < — 1, Kerr nonlinearity pro- 
vides an additional self-trapping mechanism very differ- 
ent from the above discussed two mechanisms supporting 
bell-shaped (I7I < 1) and two-peaked (7 > 1) solitons. It 
is worth stressing again that, since the two saddle points 
occur where the effective nonlinear permittivity vanishes 
^ £ ( NL ) _ Q) ) the overall discussed flat-top soliton phe- 
nomenology can be uniquely observed by means of the 
extreme nonlinear metamaterials we are considering in 
this paper. 



V. NONLINEAR METAMATERIAL DESIGN 
AND SOLITON FEASIBILITY 

The possibility of observing both two-peaked and flat- 
top solitons is substantially based on the availability of 
nonlinear Kerr media characterized by a very small di- 
electric permittivity and a parameter 7 in the ranges 7 < 
— 1 and 7 > 1, respectively. In order to discuss the fea- 
sibility of metamaterials exhibiting such unconventional 
properties, consider the composite structure reported in 
Fig. 5 and consisting of alternating, along the y— axis, 
two kind of nonmagnetic layers of thicknesses d\ and c?2, 
respectively. We here consider monochromatic radiation 
whose free-space wavelength is 1550 nm and whose elec- 
tromagnetic TM field geometry is also reported in Fig. 5 
from which it is evident that such an electromagnetic field 
experiences an isotropic response. The layers of kind 1 
are filled with a suitable semiconductor which can be 
doped to lower its plasma frequency to the point that, 
at the considered wavelength, its dielectric permittivity 
has a negative real part and a considerably small imagi- 
nary part [lH (Zinc oxide is a good candidate that can be 
doped sufficiently high, for example with aluminium or 
gallium, to satisfy these requirements [13); besides, semi- 
conductors are characterized by a standard nonresonant 
electronic Kerr response (for which 7 = 0.5). Therefore, 
for the layers of kind 1, we can choose t\ = —8 + O.li, 
Xi = -5 • 10- 20 m 2 /V 2 and 71 = 0.5 [U EH- The lay- 
ers of kind 2 have to be filled with a nonlinear medium 
whose Kerr mechanism is different from the nonresonant 
electronic one (otherwise the whole structure would ex- 
hibit 7 = 0.5) and they have to host a gain mechanism 
for compensating losses produced by the semiconductor 
(negative dielectric) layers of kind 1. Therefore, for the 
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FIG. 5: (Color online) Composite nonlinear metamaterial for 
achieving e C 1 and 7 > 1 or 7 < — 1. The layers are staked 
along the magnetic field direction of the impinging TM wave. 
Layers of kind 1 are filled with a doped semiconductor (as 
ZnO:Al or ZnO:Ga) whose doping is tuned to largely reduce 
the plasma frequency and hence to provide a negative dielec- 
tric with small absorption. Layers of kind 2 are filled with 
a nonlinear medium characterized by molecular orientation 
Kerr response and hosting dispersed semiconductor quantum 
dots. 



layers of kind 2, we choose a background medium for 
which e 6 = 2.66, Xb = 3 • 10- 20 m 2 /V 2 and j b = 3 (the 
nonlinear parameters being typical of media character- 
ized by molecular orientation Kerr response [HI), hosting 
semiconductor nanoparticles (optically pumped quantum 
dots) for which e QD = 11.8 - b.7i, xqd = 10~ 20 m 2 /V 2 
and jqd = 0.5 (these parameters characterizing suitable 
InAs/GaAs quantum dots as discussed in Ref7[22l|). In 
order to obtain the effective electromagnetic linear and 
nonlinear parameters of the whole structure of Fig. 5, we 
exploit a two-step homogenization approach |23j ] where, 
first, the effective parameters of layers of kind 2 are eval- 
uated from the background medium and quantum dots 
properties and, second, the effective parameters of the 
whole structure are obtained from those of the layers of 
kinds 1 and 2. Layers of kind 2 can be homogenized 
by exploiting standard Maxwell Garnett approach (for 
evaluating the effective linear permittivity) and its non- 
linear extension proposed by J.E. Sipe and R. Boyd in 
Ref.pjJ (for deducing the effective nonlinear parameters 
of the composite) so that, choosing the macroscopic vol- 
ume filling fraction of the quantum dots / = 0.03, we 
obtain for the layers of kind 2, e 2 = 2.79 — 0.03i, X2 = 
(3.61 • 10- 20 - 1.6 • lCT 21 i) m 2 /V 2 and 72 = 2.68+0.07?. 
If the period d = d% + d,2 (see Fig. 5) is much smaller than 
the wavelength, effective linear [lij and nonlinear [25[ 
parameters can be introduced to describe the overall lay- 
ered composite electromagnetic response so that, choos- 
ing d\ = 52 nm and c?2 = 148 nm, we obtain e — 1.5 • 
10~ 3 + 10~ 4 i, x = (1-38 • 10~ 20 - 1.21 • 10- 21 i) m 2 /V 2 
and 7 = 4.70 -I- 0.32z (where the nonlinear parameters 
have been evaluated exploiting the approach of Ref.|j|). 
Note that Re(e) = 0.0015 < 1, Re(j) = 4.70 > 3 (dif- 
ferent both from the standard values 0, 0.5, 3 and those 
achieved in Ref.JlT}) and, remarkably, that all the imagi- 



nary parts can be neglected so that the conditions for ob- 
serving the above discussed two-peaked solitons are met. 
It is worth stressing that the fulfilling of all these condi- 
tions has been achieved by exploiting both the number 
of available constituent materials and, most importantly, 
the freedom of tuning the parameters characterizing the 
composite (semiconductor doping, nanoparticle size and 
volume filling fraction, layers' thicknesses, etc.). A two- 
peaked soliton, propagating through the above consid- 
ered medium and characterized by u 2 — 0.5, requires 
the realistic and feasible peak intensity 775 MW/cm 2 to 
be excited, value which should be compared with the 
unphysical peak intensity 1.93 • 10 7 MW/cm 2 required 
for exciting a similar (i.e. with u 2 = 0.5) two-peaked 
soliton trough carbon disulfide (for which ecs 2 = 2.65, 
Xc s 2 = 3 • 10- 20 m 2 /V 2 and 7cS2 = 3). Therefore, the 
extremely high nonlinear behavior of the considered two- 
peaked soliton would require, in conventional media, op- 
tical intensities so high to prevent their observation both 
because the nonlinear response would depart from the 
simple Kerr model and, most importantly, because the 
very high intensity would literally damage the sample. 

If in the above discussed composite structure the non- 
linear susceptibility of layers of kind 1 is changed into 
Xi = — 2 • 10~ 19 to 2 /V 2 , by following the above de- 
scribed two-step homogenization approach, we obtain 
for the overall composite medium e = 1.5 • 10~ 3 + 
1 . 10" 4 i, x = (-2.5 • 10~ 20 - 1.21 ■ 10~ 21 i) m 2 /V 2 and 
7 = — 1.83 + 0.14i. Therefore the dielectric permit- 
tivity is very small, the overall nonlinear character is 
defocusing (since Re{x) < 0) and, most importantly, 
i?e(7) < — 1 so that all the conditions required for observ- 
ing flat-top solitons are met. A flat-top soliton propagat- 
ing through the considered medium and characterized by 
u 2 = 0.99a(7) = 0.32 requires the feasible optical inten- 
sity 119MW/ cm 2 to be excited. 



VI. CONCLUSION 

In conclusion we have shown that a medium character- 
ized by the constituent relations of Eqs.([2]), with |e| <C 1 
and 7 assuming values different from the three available 
in nature (i.e. 0, 0.5 and 3), is able to support the propa- 
gation of two-peaked (for 7 > 1) and flat-top (for 7 < — 1) 
solitons. These two families of electromagnetic (perfect) 
solitons are supported by two distinct self-trapping mech- 
anisms whose main difference from the standard self- 
focusing process lies in the availability of critical points 
where the effective nonlinear response vanishes (condi- 
tion made possible by the very small value of the lin- 
ear permittivity), points able to deeply affect the soliton 
shape. Since the whole predicted soliton phenomenology 
is based on unconventional values of 7 we have proposed 
a realistic and feasible composite metamaterial exhibit- 
ing such exotic nonlinear properties. It is worth stressing 
that, as a consequence of the small value of the linear di- 
electric permittivity, the intensities required for exciting 
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the considered solitons (characterized by highly nonlin- 
ear and highly nonparaxial behaviors) are physically ac- 
cessible. Therefore the soliton benchmark considered in 
the present paper allows us to argue that the two novel 
self-trapping mechanisms can be efficiently exploited for 



fine-tuning the shape of an optical beam, shape control 
which can even become active by inserting electro-optic 
or liquid-crystal based constituents within the nonlinear 
metamaterial. 
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